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ABSTRACT 

Recent observations indicate that many if not all galaxies host massive central black holes 
(BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) 
on their actions as gravitational lenses. When lenses are modelled as singular isothermal el- 
lipsoids, binary black holes change the critical curves and caustics differently as a function 
of distance. Each black hole can in principle create at least one additional image, which, if 
observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative 
distribution of magnification for images created by black holes, we find that the cross section 
for at least one such additional image to have a magnification larger than 10~ 5 is compa- 
rable to the cross section for producing multiple-images in singular isothermal lenses. Such 
additional images may be detectable with high-resolution and large dynamic range maps of 
multiply-imaged systems from future facilities, such as the SKA. The probability of detecting 
at least one image (two images) with magnification above 10~ 3 is ~ 0.2/bh (~ 0.05/bh) 
in a multiply-imaged lens system, where /bh is the fraction of galaxies housing binary black 
holes. We also study the effects of SMBBHs on the core images when galaxies have shal- 
lower central density profiles (modelled as non-singular isothermal ellipsoids). We find that 
the cross section of the usually faint core images is further suppressed by SMBBHs. Thus 
their presence should also be taken into account when one constrains the core radius from the 
lack of central images in gravitational lenses. 



Key words: cosmology: theory 
hole physics 



galaxies: formation - gravitational lensing: strong - black 



1 INTRODUCTION 

Recent observations suggest that many if not all nearby 
galaxies host massive central black holes. Empirical correla- 
tions have been discovered between the mass of the super- 
massive black hole (SM BH) and va r ious galaxy properties 
such as the bulge mass dLaorl 1200 ll; iMarconi & Hun3 120031 ; 
lHaring & Rixl 120041: iNovak et al.1 l2006t ISokeJ |2009|). veloc- 
ity dispersion (IFerrarese & Merrittl 120001 ; iGebhardt et alJ 2000 ): 
Tremaine etal] 120021: iNipoti et alj [2003: Rob ertson et al.l I2OO6 



limi t theorem in galaxy mergers with no significant physical mean- 
ing jPendl2007l) . 

Gravitational lensing is an independent, mass-based method 
to probe SMBHs. The lensing effect of single central and off- 
centre SMBHs have been studied previously (Mao et al.1 uOOll; 



Chenf2 003a.b; B owman et al.l2004l : lRusin et alj2005l : lMao & Wit 



Ql|; 
Vit3 



(Mago rrian et al.1 11998: McLure & Dunlo 



Hul 120081: iGraham et al.l 1201 ll; iFeoli et al.1 1201 ll). luminosity 



Graham 



20071) and concentratiortlGraham et alj 1200 ll; IGraham & Privet 



20071) . These correlations suggest that the growth of black hole is 



related to galaxy formation (e.g 



Kauffmann & Haehnel 



Monaco et al1l200Cl; IWvithe & Loeb 2002; lYu & Tremaim 



hneltl 
laine 



20021 ; Volonteri et al. 2003; Di Matte o et all 120031: iHaiman et al.1 
2004 ; (Ferrarese et aTfl2006l : lYu & Lull2008l ; (Bandara et al.ll2009l) . 



although such a relation can also be a consequence of the central 
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201 ll) . The effect of lensing by binary black holes has yet to be 
explored. Such systems are generated by the mergin g of galaxies 
ffi]|2002l ; lB~erczik et alfeOOolljohansson et alj|2009l) . Mergers be- 
tween galaxies are observed and predicted in the hierarchical struc- 
ture formation theory. The formation and evolution of binary black 
holes have been exte nsively studied using analytical and numeri- 
cal m ethods (see, e.g. jMerritt & Milosavlievicl200"5l ; IColpi & Dottil 
120091 for reviews and references therein). We briefly discuss the 
processes below. 

When two galaxies merge, the orbits of their associated black 
holes will first decay through dynamical friction. The critical sep- 
aration where the binary rotation velocity equals the velocity dis- 
persion of the host galaxy is often referred to in the literature as 
the "hardening radius" of the binary (see eq.U4t. This radius plays 
an important role in binary black hole evolution. At separations 
larger than the hardening radius, the black holes are "dressed" with 
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the inner cores of the stellar bulges belonging to their original host 
galaxies. The binary separation shrinks below the hardening radius 
due to the passage of individual stars which extract angular mo- 
mentum from the binary. These stars are removed from their orbits 
after the interaction with the binary, leading to the appearance of a 
"loss cone" in phase space. If this loss cone is not refilled, the or- 
bital decay of the binary may stall, leading to a large population of 
galaxies with SMBBHs. The stalling ra diu s is typically at several pc 
to several tens of pc (see, e .g., lYul2003 and lMerritt & Milosavljevic 
l2005tlColpi & Dottill2009l for reviews). 

Some other (e.g., gas) processes need to bring the binary black 
holes closer, so that gravitational radiation can rapidly merge the 
binary black holes into a single one. The problem of whether loss 
cones are refilled fast enough is still unsolved. It is, therefore, un- 
known how many supermassive binary black holes there are in the 
universe. Thus, any probe of this population will provide additional 
constraints on the formation and evolution of supermassive binary 
black holes. 

The purpose of this paper is to study the effects of SMBBHs 
on lensing properties using simple analytical models. We do not 
consider the effects of the inner stellar cores associated with the 
black holes. In these simple models, we show that the presence of 
SMBBHs can not only disturb the critical curves of the primary 
lens galaxy but also create additional images. Many new lenses 
will be discovered with the next generation instruments, such as 
Pan-Starr£]and LSST0. Some of the more interesting cases will be 
observed with higher resolution and larger dynamical range using 
other instruments from which additional images, if detected, may 
provide direct evidence for the existence of SMBHs or SMBBHs in 
galaxies. However, these images are usually very faint and close to 
each other. The Square Kilometer Array (SKA[f|), if equipped with a 
very long baseline, which will provide high angular resolution and 
dynamic range images (maps), may provide the be st possibility. 
If len ses have shallow (non-singular) central profiles dTruiillo et al.l 
120041) , central core images are predicted, but they may be destroye d 
by the presence of SMBBHs dMao et alj|200ll ; lRusin et~ai1l2005l) . 
We also touch upon this issue. 

The outline of this paper is as follows. In section 2, we present 
the basic lensing model of a galaxy with SMBBHs, e.i., the non- 
singular isothermal ellipsoid with SMBBHs. The classification of 
SMBBHs is also discussed in section 2. In section 3, we focus 
on critical curves and caustics of a galaxy with SMBBHs. In sec- 
tion 4, we discuss the cross sections of black hole images above 
a certain magnification threshold and estimate the probability of 
these images being detectable. In section 5, we study the influ- 
ence of SMBBHs on the core images in non-singular isothermal 
lens models. Conclusions and a discussion are given in section 
6. Throughout this paper we assume a flat ACDM cosmology 
with fi m ,o = 0.3, Q.a,o — 0.7 and a Hubble constant Ho = 
100 h km s" 1 Mpc -1 , h = 0.7. 



2 LENS MODEL WITH SMBBHS 

We model the lensing galaxy as a non-singular isothermal ellipsoid 
halo plus an SMBBH. This model includes the singular isother- 
mal ellipsoid (SIE) model as a special case. The SIE model is not 



1 http://pan-starrs. ifa.hawaii.edu/public/home. html 

2 http://www.lsst.org/lsst/scibook 



only analytically tractable but also consistent with models of in- 
divi dual lenses, lens statistics, stellar dy n amics and X-ray galax- 
ies llFabbiancI Il989l: iMaoz & Rixl Il993j: iKochanekl Il99l Il996t 
Groein & Naravan 1996a, b e iRix et al.l 19971; iTreu & Koopman 



200 



Rusin et al 



2003 



. .Keetoij |2003l". 
Koopmans et al. 2006] : Gavazzi et al.l 
Czoske et al.l 120081 : iDve et alj 120081 : 



Rusin & Kochane 



2007; iParkeret 



Tu et al 



2009h . 



lck 
al. 



2005 



2007 



In other 

wor ds, the core radii ar e expected to be small in e lliptical galax- 
ies dRusin &Mal |200ll; iRusin & Tegmarkl 1200 ll; iKeetonl 120031 : 
IWinn et al.l |2004 bovce et all 120061 : IZhang et al] |2007|) . In §2.1 
we outline the lensing basics for a non-singular isothermal ellip- 
soid plus a SMBBHs. In §2.2, we focus on the classification of 
SMBBHs, which is important for understanding their influence on 
gravitational lenses. 



2.1 Non-singular Isothermal Lens Model with SMBBHs 

The dimensionless surface mass density distribution of a non- 
singular isothermal ellipsoid is given by 



l 



l 



(1) 



where r c is the core radius, q is the axis ratio, E cl - = 
c 2 D B /(4irGDdDa s ) is the critical surface density, Da, D B are an- 
gular diameter distances from the observer to the lens and source, 
respectively, and Da s is the angular diameter distance from the lens 
to the source. All the lengths (xi,X2,r c ) are expressed in units of 
the critical radius, i? cr , which is also called the Einstein radius. 



Rd = -Dd#E,sis, 



4ty 



/<Tv\ 2 -Dds 

V c ) D s 



(2) 



where #e,sis is the angle subtended by the critical radius on the sky 
(#e,sis ~ 0.2 — 3 arcsec for typical lens galaxies), and the veloc- 
ity dispersion <j v is related to, but not necessarily identical to, the 
observable line-of sight velocity dispersion. We shall ignore this 
complication in our analysis and treat it as a parameter. For pur- 
poses of illustration, the source is taken to be at redshift 2, and the 
lens is at redshift 0.5. The velocity dispersion is a v — 200 kms -1 , 
and axis ratio q — 0.7, which is th e most probable axis ratio of 
early-type galaxies dChoi et al.ll2007l) . 

The lensing properties of the isothermal ellipsoid have 
been given by seve r al authors (e.g., Kassiola & K ovnetl 1 19931 : 
iKormann etlrilll994l ; iKeeton & Kochaneklll998l) . The lens equa- 
tion including a SMBBHs is given by 



2/i = xi 



\/5 , -l / \A - Q Xl 
= tan | ■;— 

$ + r c /q 

Xl — Vl 

— m,2 



Xl - Ul 



1/2 = X2 



— mi 



^ =tanh -x ( v/W* 2 

$ + r c q 
x 2 -v 2 



xTl^q 

x 2 — u 2 



7712 



(3) 



(4) 



http://www. skatelescope.org/ 



where <E> 2 = q 2 x'{ + x\ + r 2 , and mi, m 2 are the dimension- 
less mass of the two BHs respectively. We label the two BHs of 
SMBBHs as 'a' and '6'. r a , rh are the dimensionless projected dis- 
tances from the images to the black hole 'a', 'b' on the lens plane. 
(ui,u 2 ), (vi,v 2 ) are the coordinates of 'a' and '6' on the lens 
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plane 



r a = y (xi — ui) 2 + (x2 — U2) 2 
Tb = \/ (x'l - Wip + (X2 - v 2 ) 2 . 



We also define 



m = mi + 7712, 



Mbh 



ncr v R c 
G 



(5) 
(6) 
(7) 

(8) 
(9) 



where A/bh is the total mass of the supermassive binary black holes. 
Physically, A/ cr is the mass of the galaxy contained within a cylin- 
der with radius R cr , hence m is the ratio of the total mass of the 
SMBBHs to the projected mass of the galaxy within R CT . We as- 
sume that the correlation of the total mass of SMBBHs and velocity 
dispersion is the same as that for a single black hole, which has been 
studied by many authors lIFe rrarese & Merritt 2000 ; Gebhar dt et al.l 
20001 : iTremaine et alj 120021: iNipoti et alj 120031 : iRobertson et~al] 
200alHull2008l : lGultekin et alj|2009l). In our p aper, we use the re- 
cent correlation found bv lGultekin et alj d2009h 



Af b 



1O 823 M 



ct v /200 kms" 



(10) 



for elliptical galaxies, which dominate the lensing cross-sections. 
The dimensionless total mass of the SMBBHs is thus given by 



m = 2.5 x 10- 5 h 



(11) 



,200 kms" 1 , 

Notice that m has little dependence on the velocity dispersion, al- 
though we caution that there is substantial scatter around this cor- 
relation. The magnification (/i) is given by 

-1 _ dyj_dy2_ _ dy^dy^ 
dxi 8x2 dx2 dxi 
Throughout this paper, we quote the absolute magnification only 
(without parity). 

Including the source position, the lens model has 1 1 degrees of 
freedom {y, u, v, mi/ms, <J v ,q, r c ), even if we ignore the scatter 
in the correlation between m and<r v • The parameter space is large, 
and thus in this paper we limit ourselves to illustrative examples. 
Also notice that if we set the core size r c = 0, the non-singular 
isothermal model becomes a singular isothermal ellipsoid model, 
for which the influence of SMBBHs will be studied in §3 and §4. 

The lens equation has to be solved numerically to yield the 
image positions; their magnifications can then be found through 
Eq. J 1 2b . We can, however, study the influence of SMBBHs 
on gravitational lenses using critical curves and caustics. Criti- 
cal curves are image positions whose magnifications are infinite 



0). They map into caustics in the source plane. Caustics 



mark discontinuities in the number of images, so we can study 
the differences in critical curves and caustics, between lens mod- 
els with and without SMBBHs to understand their lensing effects 
(see §3). In order to see whether additional images created by 
black holes are observable, we examine the cross section when 
these images have magnifications above a certain threshold, e.g. 
jttth > 1(T 3 , 10 -4 , and 10 -5 (see §4 and §5). 



2.2 Classification of SMBBHs 

We classify each SMBBH in terms of the separation of its two 
members. First, we calculate the condition that the rotation velocity 



of supermassive binary black holes is equal to the velocity disper- 
sion. In this case, the separation of the supermassive binary black 
holes is called the "hardening" radius (dhard), 

GM bb 



4d ha 



2at 



ard 



GmM c: 
8a 2 . 



(13) 



(14) 



where M cr and m are defined in eq. (|9). The hardening radius dhard 
can be estimated as 

/ CTv \ 1-96 
dhard « 3.53 pc — — . (15) 



If the separation is smaller than this radius, we call these SMBBHs 
'hard', otherwise 'soft'. As discussed in the introduction, this may 
also be the most proba ble separation of SMBBHs (several pc ~ 
several lOpc, lYull2002h if the "loss cone" is not refilled. 

The ratio of the "hardening" radius to the Einstein radius of its 
host galaxy is 



rfhard Km 4 / O v 

__ = _^ 6 . 87X 10 (^^31 



(16) 



Two other ratios are relevant: the Einstein radius of BH to the 
Einstein radius of its host galaxy, and the "hardening" radius to the 
Einstein radius of the BH's. Again assuming a lens redshift of 0.5 
and a source redshift of 2, the first ratio is given by 

-0.02 



771 S3 4.18 X 10" 



CB.SIS 

and the second by 



200 kms" 1 



1.64 x 10" 



(17) 



(18) 



,200 kms" 1 
Both quantities are on the order of a few percent. 



3 CRITICAL CURVES AND CAUSTICS OF SINGULAR 
ISOTHERMAL ELLIPSOID LENS WITH SMBBHS 

As discussed in the previous section, critical curves represent the 
positions of images of infinite magnification, while caustics mark 
discontinuities in the number of images. For singular potentials, 
there may also exist so-called pseudo-caustics: when a source 
crosses pseu do-caustic, the image num ber changes by one, rather 
than by two fevans & Wilkinsonll 19981) . In this section, we study 
how the SMBBHs affect the critical curves and caustics of singular 
isothermal gravitational lenses (i.e., r c = 0). 

We show the critical curves and caustics of a singular isother- 
mal ellipsoid model with central SMBBHs at different separations 
in Figs. Q] and [2] respectively. The separation (dbh) is always ex- 
pressed in units of the Einstein radius, which corresponds to about 
5.1 kpc for our example. As can be seen, there are some critical 
curves near the BHs (which for convenience we will call them 
black hole critical curves) when the separation of SMBBHs is 
< dbh < 1.6. The black hole critical curves are a single contin- 
uous curve when dbh <■ 0.05, while they become several disjoint 
curves when dbh > 0.05 (Fig. [T}. These behaviours are also re- 
flected in the caustics (Fig.|2]l- When the separation becomes even 
larger (dbh ~ 0.5), the black hole critical curves become smaller 
and smaller. However, when the separation of SMBBHs becomes 
very large (dbh ^ 1-6), the critical curve of the primary lens can 
be disturbed by BHs and the black hole critical curves are merg- 
ing into the primary critical curves, so we have peculiar caustics as 
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shown in the bottom right panel of Fig.[2j Such large separations are 
not ex pected to be common in SMBBHs (Yu 200a: IColpi&Dottil 
l2009h . 

We also show some examples of images in Fig. [JJ and the po- 
sitions of corresponding sources in Fig. [2] we can see that each 
black hole of SMBBHs can generate BH-images near itself (Fig.[JJl, 
if the source is located in special positions (Fig. [2j, e.g. near the 
pseudo-caustics, there can be multiple BH-images for each black 
hole. When the separation of the SMBBHs is zero (i.e., a single 
BH), and the source is outside the primary elliptical critical curve 
(also the pseudo-caustic, see panel (a) of Fig. [2}, there is a BH- 
image near the SMBBHs, which is shown in panel (a) of Fig.[JJ As 
is well known for a singular isothermal ellipsoid lens, if the source 
is located outside the pseudo-caustic there are no multiple-images, 
so the central image must be generated by the SMBH. When the 
separations are not large, as illustrated in panels (b) and (c) of 
Fig.[JJ the black hole critical curve is continuous, and we find that 
there are three BH images for these two cases of source positions, 
as shown in panels (b) and (c) of Fig. [2] In panel (d) of Fig. [JJ 
the SMBBH separation is larger than 0.05 and the black hole crit- 
ical curves as well as the caustics become disjoint (see panel (d) 
of Fig. O. The black hole close to the source generates one BH- 
image, while the black hole farther away from the source generates 
three BH-images. There is also another image around the position 
(x = —0.081,2/ = —0.051) associated with the SIE lens. The 
cases in panels (e), (f), and (g) of Fig. [JJ are similar: BH images 
are generated close to each BH in addition to a macro image due to 
the SIE lens close to the primary critical curve. The most unusual 
case is shown in the bottom right panels of Fig. [JJ and Fig. [2] The 
separation is such large that the two black holes distort the primary 
critical curve, 8 images are created in total, including 4 BH-images 
(see panel (h) of Fig. [TJ, when the source is located in a special 
position (see panel (h) of Fig. |2j- Some of these BH-images are 
bright (\fj,\ > 0.425) enough to be potentially detectable, but the 
probability may be low because a large separation (11.6 kpc for 
our illustrative example) is required, and most black holes may not 
be at such large distances. 

In the other extreme, we have very hard SMBBHs (i.e., with 
very small separations, d < 10~ 6 in the Einstein radius). In prin- 
ciple, the rapid rotation of the binary black hole may lead to vari- 
ations in the magnification. However, for this to be observable, its 
timescale needs to be relatively short, T/4 < 10 yr, where the pe- 
riod T = (d s /(GM hh )) 1/2 . This requires a separation d < 0.04pc 
for a total black hole mass of 1.7 x 10 s M@ corresponding to 
(j v = 200km/s. This separation is even smaller than dhard, and 
thus the binary black holes will essentially appear as a single one 
for lensing purposes. We conclude that in general the binary rota- 
tion effect will be difficult to detect using current or even future 
facilities. 



4 CROSS SECTIONS AND PROBABILITY 

In this section, we investigate BH-images with magnifications 
above several thresholds: /i t h = 10~ 3 , 10 -4 and 10 -5 ; such im- 
ages are potentially detectable (see the discussion). We calculate 
cross-sections in the singular isothermal ellipsoid model for pro- 
ducing at least one or two such BH-images. 

The cross-sections are calculated by constructing the magni- 
fication map (in the image plane) and identifying the region where 
the magnification is greater than some minimum value /i t h- We then 
map this region onto the source plane, from which we calculate the 
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Figure 3. Cross-section of BH-images of SIE-SMBBHs model as a func- 
tion of separation d^h °f SMBBHs with different thresholds of magnifica- 
tion. The red, green and blue curves correspond to fj, t ^ = 10~ 3 , 10 -4 and 
10 — 5 (from bottom to top) respectively. All cases have m = 2.5 X 10~ 3 h, 
and rai/m2 = 1. Each galaxy is modelled as a singular isothermal ellip- 
soid with velocity dispersion ct v = 200 km s _1 , and axis ratio q = 0.7. 
/^th i s me threshold of the magnification for BH images. 



source cross-section. The cross section in the case of having at least 
two BH-images above a magnification threshold reflects the over- 
lapping region where each member of the SMBBHs satisfies the 
condition. To do this, we first calculate the cross section for the 
single BH case and then for binary BHs. The overlapping cross- 
section is obtained by subtracting the cross-section of the binary 
BHs from the total cross-section of the two single black holes. 

Fig [3] shows the cross-sections as a function of the separa- 
tion between the binary black holes for three thresholds /i t h = 
10" 3 , 10" 4 and 10" 5 . As expected, the higher the threshold, the 
lower the cross section of the BH-images. The peak of the cross- 
section also moves to a larger separation when the threshold in- 
creases. 

Fig.|4]illustrates the cross-section that we can detect at least 2 
BH-images, which are generated by both members of the SMBBHs 
respectively. As shown in Fig.|4] only when the separation is small 
enough, two BH-images can be generated, and the higher the fith 
is, the lower the cross section of BH-images is, as expected. 



5 SUPPRESSION OF CORE IMAGES IN A 

NON-SINGULAR ISOTHERMAL GALAXY WITH 
SMBBHS 

In this section, we investigate how core images of the non-singular 
isothermal ellipsoid (NIE) lens model are affected by central 
SMBBHs. As is well known, non-singular isothermal model can 
generate a faint core image, and the magnification and position of 
the core images are sensitive to the core size. Observationally, very 
few central images have been observed (for more, see the Discus- 
sion), wlu£h_£anbe_used to put an upper limit on the core radius 
(e.g.. iRusin & Mai 1200 ll: iRusin & Tegmarkl |200ll: iKeetonl 120031: 



IWinn et alj|2004l : iBovce et alj|2006l : IZhang et al.ll2007l) . Not sur- 
prisingly, as in the case of a single central black hole, the presence 
of SMBBHs can also demagnify and suppress the observability of 
central core images. 

We again set the velocity dispersion a v equal to 200 km s _1 , 
the axis ratio q equal to 0.7, and adopt a core size r c = 0.05. Fig. [5] 
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Figure 1. Critical curves and examples of images for the case of a SMBBH in a singular 
isothermal ellipsoid. The mass ratio for the binary system is ra\ / m,2 = 1, and the separation 
dbh is indicated in each panel, d^h is m units of rE,siS, the Einstein radius of the singular 
isothermal sphere model. The black curves show the critical curves, red triangles show the 
image positions. We also plot the positions of the SMBBHs as black points. Each case has 
m = 2.5 X 10 — 3 h, and the galaxy is modelled as a singular isothermal ellipsoid with velocity 
dispersion <r v = 200 km s _1 and axis ratio q = 0.7. The black hole image with the highest 
magnification is indicated as a larger triangle for each case. 
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Figure 2. Caustics for the case of a SMBBH in a singular isothermal galaxy for the same 
parameters as in Fig.[T] Some example source positions are shown with corresponding images 
shown in Fig. [T] The green curves show the caustics, while the blue triangles show positions 
of the sources. We also plot the positions of the SMBBHs in these panels (black points). All 
SMBBHs cases have m = 2.5 X 10 — 3 h, and the galaxy is modelled as a singular isothermal 
ellipsoid with velocity dispersion <j v = 200 km s _1 and axis ratio q = 0.7. 
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Figure 4. Similar to Fig. [3] except that cross section is for the case when 
both BHs create BH-images with magnification greater than three thresh- 
olds as in Fig.rj All SMBBHs cases have m = 2.5 X 10 3 h, and the mass 
ratio m\/m% = 1. Each galaxy is modelled as singular isothermal ellip- 
soids with velocity dispersion <r v = 200 kms" 1 and axis ratio q = 0.7. 
/^th is tne threshold of the magnification. 



shows the cumulative distribution function for the magnification of 
core images (jUcore). SMBBHs suppress the faint end of the distri- 
bution, leaving the bright end largely unaffected. A smaller separa- 
tion will suppress the faint end of the distribution more effectively 
than a larger one. For dbh i$ 0.04, the suppression is on the or- 
der of 12%, and for dbh = 0.20 the suppression is on the order 
of 3%. Non-equal mass SMBBHs lead to smaller variations be- 
tween different separations than equal mass systems. For example, 
for dbh = 0.20 with mi/mj = 3, the suppression is on the order 
of 7%, 

Fig.[6]shows the cumulative distribution function for the mag- 
nification of core images in non-singular isothermal lens with more 
massive SMBBHs (m = 0.01). As expected, the cross section of 
core images are suppressed more at the faint end of the distribution. 
To summarize, if double black holes exist in multiply-imaged sys- 
tems, then they will lead to small changes in the constraints on the 
central mass distributions in lenses. 
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Figure 5. Cumulative distribution functions for the magnification of core- 
images {ficore) in the non-singular isothermal halo with SMBBHs. The 
heavy lines show the case without BHs. The light lines show SMBBHs 
cases with separations d hh = 0.04, 0.08, 0.20. All SMBBHs cases 
have m = 2.5 X 10~ 3 h, and m\jmi = 1 (top), mi/7712 = 3 (bot- 
tom). Each galaxy is modelled as a non-singular isothermal ellipsoid with 
core radius r c = 0.05, velocity dispersion a v = 200 kms -1 and axis 
ratio q = 0.7. 



6 DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the lensing configuration due to su- 
permassive binary black holes (SMBBHs). We show typical exam- 
ples of critical curves, caustics, and image configurations. Similar 
to a single black hole, SMBBHs can create additional images close 
to them. While we have adopted illustrative values for the axis ratio 
and orientation, we have also explored other values and found no 
significant dependence on these parameters. 

For BH-images to be observable, they have to be bright 
enough to be detected. We have examined the cross-sections for 
producing BH-images with magnification above the thresholds, 
fXth = 10 -3 , 10 -4 and 10 -5 relative to that of producing multiple- 
images, 7r(9 ESIS . We write this ratio as TZbh- The probability of be- 
ing able to detect BH-images in multiply-imaged systems is given 
by Pbh(> m) ~ Kbh(> h) /bh, where /bh is the fraction of 
galaxies with SMBBHs. 

The values of 7?-bh can be read off from Fig. [3] Thus, the ob- 
servational probability of the BH-images in a multiply-imaged sys- 
tem for fi th = 10~ 3 , 10~ 4 and 10~ 5 , are about 0.2/bh, 0.6/bh 



and 1.4/bh respectively. The probabilities, where both BHs gen- 
erate at least one BH-image for the same thresholds, are 0.05/bh, 
0.15/bh and 0.3/bh respectively, i.e., approximately a factor of 
four lower. 

As shown in Fig.Q] BH-images are very close to the BHs, so 
the separation between the brighter BH-images is approximately 
the separation of the SMBBH, of which the most probable value is 
~ 10 -4 in units of .Re, sis (see eg. 116b. Thus the resolution must 
be better than « 10 -4 arcsec, which is already achievable by VLBI 
techniques. On the other hand, we need very large dynamic range 
to detect the BH-images. When the magnification of BH-image is 
10 -3 , the dynamic range required is /imax/10 -3 , where fim^ is 
the magnification of the brightest image. For most cases, the largest 
magnification /x max is about a few to ten. Thus if we want to de- 
tect a BH-image whose magnification is 10 -3 , as a conservative 
estimate, the dynamic range needs to be > 10 4 . 

It is interesting to speculate what we can learn if we do ob- 
serve two BH-images. We have 8 constraints (from the positions 
of two black hole images and two macro-images produced by the 
singular isothermal ellipsoid lens) while we have at least 10 param- 
eters even if the core size is taken to be zero. Clearly the system 



8 Li, Mao, Gao, Loeb & Di Stefano 




-5 -4 -3 -2 -1 



Magnification of Core Image log 10 (n core ) 




Magnification of Core Image log 10 (n core ) 

Figure 6. Similar to Fig. [3] except that the total dimensionless mass of the 
SMBBHs is m = 0.01. 



is under-constrained. The number of parameters will be even larger 
in more complicated models, so we will not be able to determine 
the parameters uniquely. Furthermore, multiple BH images can als o 
be produced by a single off-centre black hole dMao&Witfeoill) . 
which may complicate the interpretation. 

The presence of SMBBHs can suppress the faint end of the 
cumulative distribution for the magnification of core images, while 
leaving the bright end largely unaffected. Their effects will need to 
be accounted for in the constraint on the central mass profiles (e.g., 
core radius). There are p resently two know n lenses with a core im- 
age (PMN J1632 -0033, lwinn et al.l2004l and SDSS J1004+41 12, 
llnada et al.ll2008h . It is not surprising that SDSS J1004+4112 has 
a core image, becau se it is a cluster lens with a shallow NFW 
dNavarro et aHll997T) profile. For PMN J1632-0033, t he evidence 
for th e core image has been discussed in detail by IWinn et all 
d2004h . At frequencies higher than 1.7 GHZ, the logarithmic slopes 
of flux density ratio vs. frequency for the three images are en- 
tirely consistent with the third image (C) being the elusive and long 
sought-after central image. At 1.7 GHZ, image C is fainter than ex- 
pected, which may be due to absorption and scintillation through 
the dense lens galaxy centre. In addition, they predict a fourth im- 
age induced by the central BH at < 10% level of the flux of image 
C, which, if detected wit h VLBI technique s, will provide a mea- 
surement of the BH mass dWinn et ai1l2004l) . 

To summarize, gravitational lensing can in principle be used 
to detect single BHs and SMBBHs in galaxies through the extra im- 
ages they create. However, these images are usually very faint and 



close to each other, so they pose challenges for current instruments 
both in terms of resolution and sensitivity. For example, VLBI 
techniques may have sufficient resolution, but the dynamical range 
achievable currently may be insufficient to detect multiple, faint BH 
images. Another complication may arise because of the confusion 
of central images with radio emission from the lens galaxy. How- 
ever, most lensing galaxies are ellipticals and so their central AGNs 
may be weak. Even so, at very high sensitivity/dynamical range, 
co nfusion with centr al A GNs may still be a n issue. As discussed 
by I Winn et al. | d2004l) and lRusinetall d2005t) . we can use the usual 
tests - a common spectrum (or flux density ratio vs. frequency), 
surface brightness and correlated variability (time delays) - to dif- 
ferentiate central images from an AGNs in the l ens galaxy. For ex- 
ample, for B2108+213 dMcKean et aT]|2005l) . iMore et al.l d2008h 
compared the spectrum of the central radio source with those of 
lensed images and concluded that the central source is an AGN 
rather than a lensed image. Complications due to absorption and 
scintillation can be overcome by observing at high frequencies in 
the radio (since their effects scale as v~ 2 ). Furthermore, if multiple 
BH images are discovered, then the confusion may be less of an 
issue since there is likely only one central AGN. 

Future surveys using optical telescopes, such as Pan-Starrs and 
LSST, can provide a much larger sample of lenses, which may be 
used to identify particularly promising cases for further studies. Fu- 
ture generation of instruments, in particular the Square Kilometer 
Array, will provide very high-contrast (> 10 6 ) and high-resolution 
(< 10 -4 arcsec) imaging capabilities. We remain cautiously op- 
timistic that binary black holes can be independently discovered 
through careful observations of multiply-imaged systems, espe- 
cially in the radio. 
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